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The magnetization of a two-dimensional ferromagnetic Heisenberg model, which represents a 
quantum Hall system at filling factor v = 1, is calculated employing a large TV Schwinger boson 
ON \ approach. Corrections of order 1/N to the mean field (N = oo) results for both the SU(iV) and 

the O(N) generalization of the bosonized model are presented. The calculations are discussed in 
detail and the results are compared with quantum Monte Carlo simulations as well as with recent 
experiments. The SU(iV) model describes both Monte Carlo and experimental data well at low 
temperatures, whereas the O(N) model is much better at moderate and high temperatures. 

o 

(N| ; PACS numbers: 73.40.Hm, 75.10.Jm, 75.30.Ds 

■ T. INTRODUCTION 

^ ■ Progress in materials synthesis has allowed to study a variety of two-dimensional (2D) systems such as thin films, 
surfaces, and semiconductor quantum wells. These systems as well as the nearly 2D cuprates have lead to much 
interest in 2D quantum magnetism. It has-been found that 2D electron gases in quantum wells in the quantum Hall 
regime are novel itinerant ferromagnets iSu The strong external magnetic field quenches the kinetic energy, leading 
to widely separated Landau levels, but because of band structure effects it couples only weakly to the electron spins. 
Thus low-energy spin fluctuations play an important role. I-,,-, 

These 2D continuum ferromagnets exhibit topological excitations called skyrmionaxcl in analogy to the Skyrme 
model of nuclear physics Em In the quantum Hall system these excitations carry electrical chargeffl At filling factor 
v = 1, i.e., if the spin-up states in the lowest Landau level are just filled, skyrmions only appear as thermal excitations 
i of the form otskyrmion-antiskyrmion pairs. At filling factors away from unity, however, skyrmions appear even in the 

ground state.0 At all filling factors, low-energy spin fluctuations are also present—The combination of spin fluctuations 
^ — 1 ■ and skyrmions dramatically alters the magnetizatiomBt] and the specific heatJld 

For a quantitative understanding it is useful to first study the case of v — 1 to isolate the effect of low-energy spin 
fluctuations, which are expected to be well described by a Heisenberg model, at least at low enough temperatures. 
At higher temperatures higher oijder gradient terms neglected in the Heisenberg model could become important. 
Renormalization group arguments^ show that in D = 2 — e dimensions the magnetization M of the quantum Hall 

■ ferromagnet at v = 1 is a universal function, M/S — f(JS 2 /T, B/T), where S is the spin, J is the exchange .coupling, 
and B is the external magnetic field. For D = 2 this universality is violated by logarithmic corrections]!^ In the 

ON i Heisenbecg model the magnetization only depends on the three dimensionless quantities S, J/T, and B/T. Read and 
SachdevM have evaluated the magnetization using SU(iV) and 0(N) Schwinget boson formulations of mean field (MF) 
theory, i.e., N — * oo, for the Heisenberg model. In a recent communication!^ we have presented analytic results for 
the leading 1/N corrections to the magnetization and results of extensive quantum Monte Carlo simulations. I»,thc 
present paper we present details of the 1/N theory. Details of the Monte Carlo simulations are given elsewhere.E^I An 
alternative micr.Q§copic approach which includes spin-wave corrections to the electronic self-energy has also recently 
been developed.t-3 

Schwinger boson theoriesli3E3 have proved useful in finding MF theories, which respect the symmetry of the Hamil- 
tonian. Formal results to any order in 1/N have also been cJblainedOO However, numerically evaluating the first 
order (1/N) corrections is not an easy task. Trumpcr et al£3 have evaluated various ground-state quantities of a 
frustrated antiferromagnet in the absence of external fields. Although they are not using the large N formalism, their 
method is equivalent to a 1/N expansion to first order. 

There are a number of subtle pitfalls in the 1/N calculations, e.g., regarding normal ordering of operators. It seems 
worthwhile to present the calculations in some detail for the benefit of readers interested in using 1/N expansion 
methods. We also hope to make the physical interpretation of these theories clearer and shed some light on the level 
of accuracy of 1/7V expansions. 

In the following we give an overview of this work. First, the Heisenberg Hamiltonian is mapped ontp-.an equivalent 
boson system. There are several ways of doing this. One is the Holstein-Primakoff representation,^ which has a 
number of disadvantages, e.g., Ake square root of operators it introduces, and we do not employ it here. Instead 



we introduce Schwinger boson in two different ways. The first, presented in Sec. [II A|, makes use of the SU(2) 
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symmetry in spin space of the Heisenberg model (which is explicitly broken by an external field) . The second utilize s 
the local equivalence between the groups SU(2) and 0(3) to write down an equivalent 0(3) boson model (Sec. Iff A| ) . 
Subsequently, the two models are generalized to SU(iV) and 0(N), respectively, which contain N bosons at each site. 
At this point a remark may be in order on what we do not mean by the 0(N) model, ft is not an TV component 
vector model, e.g., an N component quantum non-linear sigma model. Rather, the spin operators are generators of 
the Lie group O(iV). Only for N = 3 are the generators antisymmetric 3x3 matrices, which are dual to (axial) 
vectors. Thus our results are not easily compared to expansions in the number of components of the spin vectors, as 
developed by GaraninEiJ for a classical system. _ . 

It is now possible to expand in 1/N as a small parameter. MF theory becomes exact for both SU(oo) and O(oo).tll 
The 1 /N expansion is a saddle-point expansion around this MF solution, not a perturbative expansion in the inter- 
action. For this reason it is, in principle, equally valid at all temperatures. Also, it does respect the symmetry of the 
Heisenberg model. This property makes even the MF results qualitatively correct. In particular, the absence of long 
range order if no external field is present is correctly predicted. After rederiving the MF magnetization in Sec. II, 



for the SU(iV) mode l and in Sec. IIIB for 0(N), we calculate the 1/N corrections using a diagrammatic approach! 
(Sees. |ll C| and [II C ). These corrections take fluctuations around the MF result into account. We will make use of 
gauge invariance to simplify our task. Here we also have to discuss the effect of normal ordering. In principle, terms 
to any order in 1/7V can be obtained in the same way. 

The system without exchange interaction can be solved exactly for any value of the spin S and for any N in both 
the SU(7V) and the 0(N) model. It can be used to check the I/TV expansion. However, the interaction introduces a 
number of additional complications. 



II. SU(iV) MODEL 
A. General considerations 



We start from a Heisenberg model with nearest-neighbor interaction on a square lattice in a constant magnetic 
field, 



// = -J^S( i ).S( J )- J B^5 z ( ? ), 



(i) 



where the sum over (ij) is over all nearest-neighbor bonds. A factor of g/iB has been absorbed into the field B. The 
total spin at each site is S; S(i) ■ S(i) = S(S + I). We express the spins in terms, Hose operators using a Schwinger 
boson representation, where two Bose fields a and b are introduced according tc 



S+ = a)b, S- = tfa, S z = (a f a - /2. 



(2) 



To restrict the Hilbert space to the physical states, the constraint a^a + b'b = 2S is introduced, which corresponds to 
S • S = S(S + 1) for the original Hamiltonian. The boson Hamiltonian is 



H=-~J2 [a\i)a(i)a^U)aU) + a\i)b(^ ti)a(j) 

+ b\i)a{i)J{j)b{3) + b\i)b(i)tf(j)b(j) 
-f 52[aHi)a(i)-tf(i)b(i) 



(3) 



neglecting a constant. For this Hamiltonian to be equivalent to the Heisenberg model, the spin operators expressed 
in terms of bosons have to have the correct commutation relations. This is easily shown to be the case. 

Utilizing the SU(2) symmetry group of the spins we write the Hamiltonian in a more compact form by first defining 
a SU(2) spin matrix 



a^a a^b 



(4) 



with the constraint Tr S = 2S. The Hamiltonian is 
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(5) 



where a z is a Pauli matrix and summation over repeated indices is implied. Here, the spin matrices S(i) should be 
infinitesimal generators of the SU(2) group, i.e., elements of the corresponding algebra. This is not the case since 
the generators are traceless. However, if we had defined S as an element of the algebra, the Hamiltonian would only 
change by a constant and we use the more convenient definition (Q). 

The group SU(2) is generalized to SU(iV) for any even N. The generalization of the Hamiltonian is 



H 



N 



(6) 



(ij) 



where S and h are N x N hermitian matrices and S is subject to the constraint Tr S = NS. We choose hg — 5 a p{— l) a+1 
so that we regain the SU(2) model for N = 2. The Schwinger boson representation now requires N boson species 
b\jbp, and the constraint is 



u a ,\ — I Op 



blb n = NS. 



(7) 



We now go over to the continuum for mathematical convenience. The continuum model may actually give a better 
description of itinerant magnets but is harder to compare to Monte Carlo simulations on a lattice. Up to a constant 
we obtain 



H 



d 2 r 



J 

2N 



(d 3 S a )(d 3 Si) 



2a? a 



(8) 



where b a (r) is a continuous Bose field with the commutator [b a (r), bg(r')] = a 2 5 a p5(r — r'), dj is the two-component 
gradient, a is the lattice constant, and summation over j is implied. After bosonization we find 



H 



d 2 r 



JSidjbbidjbJ - jfbUdjbDbeidibj 



(9) 



which is normal ordered, necessary for the functional integral. This is basically the Hamiltonian of the complex 
projective CP 1 ^^ 1 model.ES We have used the fact that the lattice Hamiltonian (^) can be normal ordered trivially 
since spins at different sites commute so that Sp(i)S (j) = : Sp{i)S^{j) : , where : : denotes normal ordering. 

Now we write down the partition function as a coherent state functional integral, where the Bose fields are replaced 
by complex fields, 



Z = / D'b a D\expl -- 



n/3 
dr 



d 2 r 



C[b;X] 



(10) 



where the functional integral takes each b a (r, t) over the whole complex plane and each A(r, r) parallel to the imaginary 
axis (a constant real part is irrelevant). Here and in the following we neglect constant factors in Z. r is the imaginary 
time, (3 is the inverse temperature, and C is the Lagrangian 



C 



^b* a d Q b a + JS{d 3 b* a ){d 3 b a ) - ± b * a (dS)b fJ (d 3 b a ) 



B 
27i 



n 2 ^/3^/3^ Q 



Xb* a b a - NSX. 



(11) 



The first term is the usual Berry phase (do is the time derivative) and the last two terms come from the constraint 
using the identity 2ir5{cj)) = dxe^. A is a Lagrange multiplier at each point (r, r) . 
To decouple the quartic term we introduce a Hubbard-Stratonovich field Q(r,r): Since 



DQj cxp — 



Tui 

dr 



d 2 



r^[-iNQ 3 -{d 3 bl)b 



x[iNQ 3 -b* p (d 3 b l3 )] 



(12) 
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is independent of b a , we can multiply the partition function with this expression. Qj can be chosen real since an 
imaginary part of Qj would not couple to the b a fields because b%djbp is purely imaginary. We get 



J D 2 b a DXDQ je xp(~^ ^dT J d 2 r£[b;\,Q] 



with 



£ = ^b* a d b a + JS(d 3 b* a )(d 3 b a ) + NJQjQj 



(13) 



+ UQ j b* a (d j b a )-iJQ j (db* a )b a 
We see that Q is a gauge field: If we multiply all b a by a local phase factor, b a (r, r) — * e %e ^^b a {r, t), we reobtain 



^h a b* b a + Xb* a b a ~ NSX. (14) 



the Lagrangian (14) by letting Qj — > Qj + S djO. We know from gauge theory that Q contains more information than 



is physically relevant; we have the freedom to choose a gauge. We use a transverse gauge, 

djQj = 0. (15) 

Of course we obtain the same results if we do not fix the gauge. The gauge freedom then leads to the appearance of 
zero modes, which turn out not to enter in the magnetization. 

B. Mean field theory 

Up to this point the treatment has been exact. In the following we derive mean field (MF) results, which are 
exact for N — > oo and approximate for finite N. This approximation is not the same as standard MF theory for the 
Heisenberg model. As we will see, SU(7V) MF theory captures the low-energy spin-wave physics of the Heisenberg 
model and correctly predicts the absence of long range order at finite temperatures. 

The MF approximation is the leading order of a stationary phase approximation for the SU(iV) partition function. 
The MF solution is assumed to-be homogeneous and static, i.e., Q and A_are assumed to be constant. This assumption 
is not justified for all systems,E3 it should hold in ferromagnets, though. O The MF values Q and A are chosen in such 
a way that the MF free energy Fo has a saddle point. If we set A to its MF value the constraint (Q) is no longer 
satisfied locally but only on average. In order to diagonalize the action we introduce Fourier transforms of the b a 
fields, 

b a (r,r) = ^ J ^^e'^Xfk,^), (16) 

where iu> n = i2irn/h(3 are bosonic Matsubara frequencies. From now on summation over indices is written out. With 
Eq. ( p^ ) and the definition hp — 8 a p(— l) a+1 the MF partition function is 

Z = J D 2 b a (k,iuj n ) exp ( - NN/3JQ ■ Qa 2 

+ JVa 2 NS(3\- I 'ft^4>] , (17) 

iuj„ ) 

where M is the total number of sites and 

= (3a 2 J2 ( - ih ^n + JSk 2 a 2 - 2 JQ • k, 



a 2 



2 

We introduce a number of new symbols, 



B h a a + a 2 \) b* a (k,iu; n )b a (k,iu; n ). (18) 
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A = a z (3\ - --Q ■ Qa z , J = f3J, B='—. (19) 
o it 



Evaluation of the Gaussian integrals yields 

z Q oc e^nn n ( ^« + j5fcv 



-2JQ-ka 2 - Bh^ + A+^Q-Qa 2 J . (20) 

Writing the product as the exponential of a sum, replacing the k sum by an integral, ~~ * (Afa 2 /4ir 2 ) J d 2 k, and 
shifting k by Q/S*, we obtain 



Z oc exp (mNSA - ^ J d2k /2z2 



xln -i(3hcj n + JSk 2 a 2 - Bh a a + A ^j. (21) 

The MF partition function and thus all MF quantities only depend on A and Q through A. The saddle-point equation 
for A is d\n Zo/dA — 0, resulting in 

iuj n ex. 

x _ (22) 

-i(5huj n + JSk 2 a 2 — Bh% + A 

The Matsubara sum in this expression is not well-defined since the summands do not fall off fast enough. In writing it 
as a contour integral the contribution from closing the contour does not vanish. The usual procedure is to introduce 
a convergence factor e ±l?,/3?iaJ " and let 77 — s- afterwards. The result is ambiguous, depending on the sign in the 
exponential. Here, Eq. (E2) only has solutions for positive sign. Consequently, 



0=AAW5-^-^ J d 2 kn B {JS k 2 a 2 - Bh^+A) 



■ NNS + ^^Y\n(\-e- A +zK). (23) 



Here, ^(e) = l/(e e — 1) is the Bose function. Eventually we fincfcJ 

1 



S = — 



8irJS 



ln(l-e- A+B )+hi(l-e- A - B ) . (24) 



Equatio n (B4| ) for A can be evaluated analytically. For given A we have the freedom to choose Q, and A is then fixed 
by Eq. (|19| ). This is a consequence of gauge invariance since Eq. ( |l5| ) specifies the gauge only up to a constant. We 
choose Q = 0. (The square lattice model without continuum approximation runs into problems at this point since 
the quantity corresponding to Q shows a spurious first-order transition at the MF level.) _ 
The MF magnetization normalized so that Mq(T = 0) — S can be obtained from Eq. (pi|)liiJ 

Mr, = — z 111 Zn 

AfN{3 dB 



1 



8irJS 



ln(l - e- A+B ) - ln(l - e~ A - b ) . (25) 



Some notes are in order: (i) Equation ( pi[ ) states that the total number of "up" and "down" bosons (with h™ = 1 
and —1, respectively) is conserved, whereas Eq. (|2^) states that the magnetization is basically the difference of the 
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number of "up" and "down" bosons, (ii) The dependence of Z$ on the field B through A is irrelevant at the MF level 
since <91n.Zo/<9A = by definition. This is not the case at the 1/N level. (iiiLThe normalized magnetization Mo/S 
exhibits the universality mentioned in Sec. |: It only depends on JS 2 and Bx3 

Finally we compare the MF magnetization with the original Heisenberg model. From Eqs. ( p4| ) and ( p5| ) we obtain 
at low temperatures 

Mo - S= ln(l - e _/3B ) (26) 



up to exponentially small corrections to the field B of order of A — B = e ~ 87rJS /(l — e~P B ). However, Eq. (|2g) is just 
the magnetization of the Heisenberg model neglecting magnon interactions. This means that the SU(7V) MF theory 
captures the correct low-energy spin-wave physics. Consequently, we expect higher order corrections to be small for 
low T. 

C. 1/N corrections 

To take fluctuations in the auxiliary fields A and Q into account, we write 

A(r,r) = A + iAA(r,r), (27) 
Qi(r,T) =0 + AQ,-(r,r). (28) 

The fluctuations in A are imaginary since A has to be integrated along the imaginary axis in Eq. ([T^) . The fluctuations 
in Qj are real. They are subject to the gauge constraint in Eq. (|l5|). 

We follow the procedure outlined by AuerbachJlH The exact partition function is 

z = J DAXDAQj exp(-NS), (29) 

where the action S is expanded in a series for small fluctuations rg with rg standing for any mode AA(r, r) or AQj (r, r), 

5 =E^:v^--^> (3°) 



nl 

Tl = 



where summation over repeated field indices-Iji is here and in the following implied. On the other hand, the action 
can be written as S — So + Sdi r + Si 00 p witbllZI 

So = iTr In G -\ (31) 
1 f hf3 f 

Sdix = -T7T / dr d 2 r (NJQ ■ Q - NSX), (32) 
Nh Jo J 

Soop = — Trln(l + G v t r e ) , (33) 

where the trace sums over space, time, and boson flavor, Go is the MF bosonic Green function, and vg is a vertex 
factor coupling the fluctuation rg to two bosons. 

The first term, So, has the standard form for a non-interacting system. It stems from the k integral part of the MF 
free energy; see Eq. (^l[). The Green function can be read off from the MF partition function, 

G^(k,ioj n ) = {-ihu n + JSk 2 a 2 - Bh a j2 + a 2 A) _1 . (34) 

The second term, Sd; r , comes from the constant part of the MF free energy but also contains fluctuations in the fields 
rg which do not involve bosons. The constant part conspires with So to form the MF free energy —/3Fq = NS^-°\ The 
fluctuating part contains a first order term in AA, corresponding to the coupling of A to the constant NS in Eq. (p^|), 
and a second order term in AQj from the Q • Q term. The corresponding diagrams are shown in Fig. [I]. 

The third term, Si oop , contains the contribution of fluctuations rg coupling to bosons. It is the result of a linkcd- 
cluster expansion. By expanding the logarithm we obtain the contribution from Si oop to S^ n ' , 



G 



c(«) 
>->«. 



loop N 



1 (— l) n+1 

y2Tr(G v ei ...G v e J. (35) 



The sum ^ p runs over all permutations of the n vertices. The first few terms are shown diagrammatically in 
Fig. ||. Solid lines with arrows denote MF boson Green functions Go and the dots correspond to vertex factors vg. 
The wriggly lines are external legs re. . Disconnected diagrams are taken care of by a linked cluster expansion, which 
puts the whole series into the exponential. No internal rg lines appear since as far as the action is concerned the 
r^(r, r) are external variables. 

For n > 3 Eq. (|35|) is the only contribution, whereas and contain contributions from iSdir and <Si 00 p- The 
total first order term can be shown to vanish as it should since we are expanding around a saddle-point. 

To find the vertex factors vg we write the exact partition function Z of Eq. ( [l3|) in terms of Fourier transforms, 
where the b a dependent part of the Lagrangian is 

C" = (3a 2 (-ihun + JSk 2 a 2 - ^h a a 

a ^ ~ 

4 



/3a r 



x 



(-2 JQ(q, iv n ) ■ ka 2 + a 2 A(q, iv n )) 



x6*(k,ia; n )6 Q ,(k- q,iu> n - iv n ), (36) 

The first expression in parentheses is the inverse Green function (Gq The same prefactors have to be included in 
the vertex factors, which are the coefficients of the terms reb a b a . Consequently, 

a 2 47r 2 2 2tt 2 
VAX= 2^If^ m =AA m ' (37) 
a 2 4n 2 , 9 2n 9 . , 

VAQ > = 27T JTa? ( " 2J) 3 = ~77 2Ja r (38) 

The factor Air 2 /Afa 2 in both cases stems from the integral over q. The factor of i in v&\ comes from Eq. (fj7j). 

We now consider the expectation value (b a b a ) for any a (no summation implied). From this we obtain two important 
quantities: The average number of bosons per site n = Y^aWx^a) > and the magnetization M = N^ 1 J2 a ha(b a b a ). 
Inserting a source term A£[j Q ] = ^2, a jab%b a into the Lagrangian (fl4|), where the source current j a is constant, we 
find 



{ a a) ~ M(3a 2 Z d 3a 



(39) 



Inserting the series expansion of Eq. (pOj), evaluating the derivative, and expanding the exponential of the terms 
containing , n > 3, we obtain 



x 



N f f^l 55, (n) 

Jd^daq, £ 

' \n=0 
m— \ n — 3 




S?X r ^ ) ■ ( 4 °) 



All terms are Gaussian integrals, which can be evaluated by pairwise contraction over the fields rg. Diagrammatically, 
any contraction is represented by connecting two vertices by an RPA fluctuation propagator D — (S^)^ 1 , which we 
represent by a heavy wriggly line. 

In the next step we calculate the j a derivative of S^. The derivative basically replaces Go by —(Go) 2 so that 
we may expect it to be related to The Green function in the presence of the source term is Gg(k, iu> n ) — 

(—ihuj n + JSk 2 a 2 ~Bha/2+a 2 A + a 2 j a )~ 1 so that its derivative is OGq /dj a = — Gq^Gq • The vertex factor associated 
with j a differs from v&\ only in a factor of i, vj a — 2irAT~ 1 a 2 . With Eq. (|35|) we have 



_d_ ~(n) 



Af 1 (— l) n+2 v-^ 

= 7T- - — V Tr (G v Ja G v ei . . . G Q v ln + ... + G a v ei . . . G v joi G Q v tn ) . (41) 

J"=° p„ 
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The sum contains nnl terms and not (n + 1)1 because Vj a cannot appear to the right of V£ n . The invariance of the 
trace under cyclic rotation allows us to write this expression as a sum over all (n + 1)! permutations of the, vertices 
Vj a , v^, . . . , ve n , if we introduce a correction factor for overcounting, nn\/(n + 1)! = nj (n + 1). We obtaintH 



d c(n) 



M 1 (-1)«+ 2 ^ n.. ^_^ c( „ +1) 



— Tr {G v ja G Q v tl . . . G v tn ) = — SfJZin- ( 42 ) 



j a =Q 

Equation (17.25) in Ref. [IT] differs from this result because of different definitions of vertex factors. It follows that 



<^ = +^ / DAX DAQj [E^ffi^...^] 

x E -—{- E ~i s S..eJ^ • ■ ■ r ^ ex P ( - -o- 5 iiV* r fc ' 

m=0 ' \n=3 ' / ^ ' 

In principle we can evaluate the integral for any term in this series. The contraction of two variables gives 



(43) 



- J Dn r £l r e2 exp(-^4J 2 r^) = ^(S {2 %\ 2 , (44) 

where the r# are assumed to be real. (In fact they are real only in direct space but complex in Fourier space. We can 
use the Gaussian integral for complex fields, which has an additional factor of 2, and note that AA(— q, — iv n ) and 
AA(q, iv n ) are not independent since AA(— q, — w n ) — AA*(q, iv n ) and similarly for AQ. Thus we have to restrict 
the sum over q to one half space. The factor of 1/2 obtained in this way cancels the factor 2 from the Gaussian 
integral.) The RPA propagator D is the inverse of the matrix 6>( 2 '. We obtain all terms in the expansion ( ^3| ) by 
writing down all allowed diagrams consisting of any number of boson loops with one external j a leg, represented by 
a dashed line, and an even number of internal vertices, and connecting the latter by RPA propagators in all possible 
ways. Allowed diagrams are connected and do not contain loops with only one or two internal vertices and no external 
vertex because first order terms cancel in an expansion around a saddle point and the loop with two internal legs is 
already included in the RPA propagator. 

To figure out which terms are of which order in l/N, note that the magnetization M — N^ 1 ^2 a h^blfia) contains 
an explicit factor of 1/A. Furthermore, each loop contributes a factor of N from summation over flavors (the loop 
with the external vertex does not contain a sum but has an explicit N), and each RPA propagator contributes a factor 
of l/N; see Eq. @. 

The leading term is depicted in Fig. ||. It is of order N° in the magnetization. This term reproduces the MF 
magnetization ( p5f ) . Contributions of order 1 /N in the magnetization have to contain the same number of loops and 
propagators. The only two allowed diagrams are shown in Fig. || Similarly, we could write down the diagrams to any 
order. 

The two relevant contributions are, from Fig. ^(a), 



and from Fig. ||(b), 



1 



_c(l+l) c(3) (e(2)\-l f,c(2)\-l ( 4 a\ 



T>2 contains a sum over all three possible pairings of the four n. We have utilized the symmetry of in each 
pair of indices. The evaluation of the relevant is relegated to App. [A|. With the results found there we can write 
down the two terms in (b^ a b a ) = 2?j + T>2 i 

rr a/ \ 

(47) 




MN n B {k-B)+n B {K + B)^ W 4 NN/2J + Y,a a l) ' 
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where we have introduced new symbols, 



— Oil 



^(6g +q/2 )-^(e£_ q/2 ) 
^ -if3hv n + 2JSqk ia 2 ' 

y 2fl2 ^( £ £ +q/2 )^B( £ g_ q/2 ) 

^ -if3hv n + 2JSqk ia 2 

n B ) ( £ k+q/2) ~ n B ) ( £ k-q/2) ^ 

~* —i(3Kv n + 2JS qk\a 2 

Vfc 2 V^ )(££+q/2) ~ n } 1)(e£ - q/2) (52) 
V -i(3hv n + 2JSqk ia 2 



with = JS k 2 a 2 — Bh" + A and the derivative of the Bose function. 

Before we turn to the magnetization we look at the total number of bosons per site n = Ea(^a^a)- The MF 
contribution is rio = NS since A was chosen that way. To the next order, 



1 JVn B (S-B)+n fl (S+B) ^ / EX' , 



AAiV 2 n B ( A - B) + n B (A + B)^\ ^ a$ NN/2.J + 



= NS + 0. (53) 



We thus see explicitly that the constraint (Q) is still satisfied at the 1/N level. This is a special case of Auerbach's 
general prooO that the constraint is satisfied to any order. Thus we need not "shift the saddle point," i.e., adjust A 
so that the constraint is still satisfied. 

We now calculate the 1/N contribution to the magnetization, 

M = Mq - * V ( E ° (C1 ~ f>°? + E " (C1 -" K)< 3 ) , (54) 

2 ^irA EX MN/2J+Y.A) 



where Mq is the MF magnetization (|25|) and 



n B (A- B)-n B {A + B) 

ci = = = = — . (55) 

n B (A- B) + n B {A + B) 



Expressing the momentum sum by an integral we find 



N ^ 2 J tz\ ex ^/2j+eX/ 

= Mo "I^ / A S ( AM o(q, ^») + AMi(q, «/„)) . (56) 

The first term in the parentheses is due to fluctuations A A while the second comes from AQ. 

The same result ( pq) can be found by writing down the one-loop contribution to the free energy and taking the 
derivative with respect to magnetic field. In this way we find physical interpretations for the two 1/N diagrams: 
The term coming from the explicit B dependence of the free energy corresponds to Fig. ||(a), whereas the indirect 
dependence through A(B) corresponds to Fig. |J(b). 

Although the expression ( |56| ) for the magnetization is formally correct, great care is needed in evaluating the 
frequency sum. We will first show briefly how naive evaluation leads to a spurious divergence of the integral over 
momentum q and then present the solution to this problem. The solution involves carefully taking into account 
normal ordering of boson operators. First, we derive the contributions to the magnetization for large momentum q 
of the (in Fig. ^ vertical) RPA propagator. We express the momentum sum in <7q and er" by an integral and shift the 
variable k in the two summands by k — ► k =p q/2, 
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°'- L 4tt 2 7 " \fc|a 2 / \-ij3%v n + 2JSqk ia 2 - JSq 2 a 2 -ij3hv n + 2JSqk ia 2 + JSq 2 a 2 

xn B (JSk 2 a 2 - Bh% + A), (57) 

where the upper (lower) expression in {} pertains to erg (er"), and ctq' and c"' are obtained by replacing n B by n^ B . 
Formally, we expand for small k\ since large k are exponentially suppressed. Odd powers of k\ vanish so that 

ffQ = A^ ^ / 1 1 

0,± 47r2 m CTen \ H/^n ~ ^ 9 2 a 2 ) m+1 {-if3hp n + JS q^a?)"* 1 

x J d 2 k^ k l a2 ^(2JSqk ia 2 ) m n B {JSk 2 a 2 -Bh*+A). (58) 

This is an assymptotic series and does not converge. However, for our argument we only need the first two non- 
vanishing terms, which are well-defined. 
The leading term in J2a a o rea ds 

y a % s -JL. ( ) y i n( i _ 6 -a+^S) 

^ 4irJS \-ip?W n - JSq 2 a 2 -i/3fiu n + JS q 2 a 2 J ^ 

= AfNs( ^ V (59) 

\ —i(3fiv n — JS q 2 a 2 —i(3Tiv n + JS q 2 a 2 J 



where the last step follows from the saddle-point equation (|24|). Thus the to = term is independent of magnetic 
field. Similarly we find ^2 a (ci — h™)^' = to the same order, where we have used Eq. This result is not 

surprising since X) a ( c i — ^ l a) <J d' i s > U P to a factor, the magnetic field derivative of J2 a a o- The leading non- vanishing 
(to = 2) term is 



4tt 2 \(-ij3hv n - JSq 2 a 2 f (-i/3hv n + JSq 2 a 2 ) 3 , 
x J d 2 k (2JS qha 2 ) 2 y( Cl - h a a ) n%\JS k 2 a 2 - Bh a a + A) (60) 



and consequently, for large q, 

Ecfci-W, Ja 4 / 1 1 



AM (q,w„) = 



<?o 2it2n V -ipthun - JS" 9 2 a 2 -i/3fti/„ + JS q 2 a 2 , 

x y( Cl - h%) J d 2 kk\ n% ] {JS k 2 a 2 - Bh a a + A). (61) 



Adding the two fractions under the sum and performing the Matsubara sum we find 

V ( U ) = V ~ 2JS ^ = - coth i^-. (62) 

f£ \-i{3Tiv n ~ JSq 2 a 2 -i(3n Vn + JSq 2 a 2 ) f£ {(ifiv n ) 2 + ( JS q 2 a 2 ) 2 2 

For large momentum q we thus obtain 

X>M (q,H/„) = - i ^j^E( c i - ^) ln (! - ^ Ml °) = ci - ^. (63) 

The large momentum limit of the AQ contribution is found more easily. We have to consider AMj_(q, iv n ) for large 
momenta; see Eq. (|6|). The leading order (to = 0) term of J2p a i_ 1S negligible compared to the constant added to it. 
Furthermore, the to = term in X) a ( c i — ^■a) a ± does not vanish. Thus the leading term is 



Ja 4 



2tt 2 N V -iphv n - JS q 2 a 2 -ij3Tip n + JS q 2 a 2 



X]( Cl_/l S) / d 2 kk 2 n { £{JSk 2 a 2 -Bh«+k). (64) 
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This expression is equal to Eq. (|6l]). Consequently, for large momenta the integrand of the external momentum 
integral in Eq. ( |56| ) is 

2 (ci - Mo/5) , (65) 

which is independent of momentum. This term would lead to a strong UV divergence of the q integral in Eq. 

We now discuss the cure. As mentioned above, we have to normal order the operators before we can write the 
partition function as a functional integral. Careful treatment of ordering is often essential to resolve ambiguities in 
the expectation value of operator products at equal times. Here we are interested in (b' a (r, r)6 Q (r, r)). One way of 
dealing with these ambiguities is to split the time of the operators in the Lagrangian in such a way that the creation 
operator is at an infinitesimally later time than the annihilation operator. The time-ordered product in the definition 
of the Green function then takes care of normal ordering at equal times. 

If the action contains a term like 



1 



hp 



AS = lkl dT I d2r >] c ^ T ) b a(r,T + V )b a (r,T), (66) 



where r\ > is small, then, after Fourier transformation 

f3a 6 



AS 



2-kN 



J d 2 kd 2 q Yc a (q,iv n )e l ^ 7 >bl(k,iLj n )b a (k-(i,iuj n -iv n ). (67) 



The coefficient c a obtains a phase factor e luJ ™ n , where iuu n is the frequency of the boson created at this point. We split 
the time in this way in the exact Lagrangian ( |l4| ) as well as in the source term j a b* a b a - As a consequence, phase factors 
containing the frequency of the outgoing boson appear at all vertices. The only places where they are relevant turn out 
to be in S^ and S (2+1 \ In the first term in Eq. ftA7\) we obtain a total factor of e ^-n.v^(^+^)v^n = e ^ n r, e ™„v 



from the three vertices. The factor containing iui n is irrelevant since the sum over iuj n is unambiguous anyway. 



We are left with an overall factor of e lUnri . The second term from the symmctrization in Eq. (A7) obtains a factor 
e -iv n r\ _ jf -j-jjg .j. w0 terms are added, the terms which have the denominators squared obtain a prefactor of 2ismv n r], 
which vanishes in the limit r\ — > (the denominators are already of second order in frequency so that ambiguities or 
divergences do not arise) and the remaining terms are 

5 (3) _ (2^ l0a , V V ^ ^ ) (^ +q / 2 )-^"4 1) ^- q / 2 ) ( M m 

^AA(0 ) 0),AA(- q ,- < ,„),AA( q ^„)-^3 iV ^ a 2^2^ - lf 3hv n + 2J S qk^ ( P [ ' 

^AA(0,0),AQ 2 (- q ,- w „),AQ 2 ( q ,^) " ^3^ ^ + 2 JS qk^ ^ ( ' 

5(2+!) follows as above. 

The Matsubara sum of the leading large q term is now, instead of Eq. (p 



E 



i(3Tiv n — JS q 2 a 2 ~i(3Tiv n + JS q 2 a 2 J ^^ + 

-2JSq 2 a 2 2j3T%v n 



E 



■ COS Vr, 71 H = sm V n 71 

{phv n ) 2 + {JSq 2 a^ {(3h Vn ) 2 + {JSq 2 a*y 



, JSq 2 a 2 ,. sinh JS q 2 a 2 (l/2 ~ 7]) , JSq 2 a 2 

= -coth h lim 4 — — = 1 - coth . (70) 

2 r,^o+ sinhJ5g 2 a 2 /2 2 

The sine series exactly cancels the leading term of the naive series for large q, making the remaining expression 
exponentially small. In other words, the phase factors introduced to ensure correct operator ordering just remove the 
constant (|6^) from the integrand for large q. The exponential factors are irrelevant in all other terms, which are of 
higher order in l/w n and are thus unambiguous. 

Deriving a convergence factor e ±lUnTI is a common method to resolve the ambiguity of a Matsubara sum. What 
is unusual here is that two different factors appear for the two terms. It is easy to fall into the trap of thinking 
that one should simply add the two fractions under the frequency sum in Eq. ([32|), arguing that the sum then looks 
unambiguous. This loses an essential contribution because of the two different phase factors. 
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We utilize the above result by calculating the q integrand numerically without taking normal ordering into account 
and then subtracting the constant (|6^) explicitly. The frequency sum is expressed in terms of a contour integral. 
First we analyze the analytic structure of aft etc. in the complex v plane. If we replace the sum over k by an integral 
in Eq. (^) for <7q we see that the integrand of the k\ integral has a pole and, consequently, Cq has a branch cut. 
Furthermore, it can be shown that ^ Q <Tq does not have zeros apart from the trivial case q = 0. The quantity 
5Z a (ci — ^qVo' obviously also has a branch cut along the real axis and does not have poles. Consequently, the AA 
contribution to the magnetization, AMq, has a branch cut and no poles. The contour of integration, C, is shown in 
Fig. ||. We have 



i V AM = ^-. I dvn B {(3hv) AM 
np .■ e —^ n Zni J c 

dvn B (f3hv) AM o (0) 



1 



C,- 



— e poo 



— oo J e 



dv nsiPTiv) [— iAq(v)] 



(71) 



where ^o(^) = — 2ImAAfo(f + iS) is the spectral function of AMq and C e is a positively directed circular path of 
radius e — > around the origin. Since Aq vanishes continuously at the origin we get, after moving the last term to 
the left hand side, 

1 1 f°° 

-VAM„ = -- dvn B {PKv)A {v). (72) 

W 2lT J-oo 

A similar equation holds for the AQ contribution, AM±. We do not discuss the numerical methods in any detail 
since they are standard. We only note that it is useful to expand the Bose functions n B in <7q etc. in a geometric 
series because this allows one to perform the integrals over k analytically, thereby replacing 2D integrals by numerical 
summation of well-behaved series. 

After evaluating the frequency sum, we have to subtract the constant (|65[). Numerically the correction term is 
indeed found to cancel the constant for large q. The new leading term drops off as 1/ q 2 so that the q integral diverges 
only logarithmically. We regularize the integral by restricting it to the first Brillouin zone, i.e., by a lattice cutoff. 
We use a circular Brillouin zone. The integration over the angle of q is then trivial. 

We find that fluctuations in A and Q always decrease the magnetization, as is intuitively expected. In fact the 
magnetization to order l/N can become slightly negative. Of course, the exact magnetization cannot be negative. 
Apparently the l/N expansion does not work well for SU(2). We can force the result to be positive by putting the 
fluctuations into an exponential: In writing down the functional integral we should impose the constraint that the 
total magnetization be positive. This constraint can be implemented by writing the full magnetization as M = Moe 9 
and expanding g in powers of l/N. If M = M + M 1/N + 0{1/N) 2 then 3 = 1 + M 1/N /M + 0(1/N) 2 . Both 
expansions are equally valid to order l/N. Of course, this method seems dubious if the l/N term is not small. 

In Fig. H we plot the magnetization for S = 1/2 and the fields B/J = 0.05,0.1,0.25 as a function of T/B. At 
the lowest temperatures the SU(iV) model describes the Monte Carlo results quite well, as expected since the SU(iV) 
model captures the correct low-energy physics. The 1/iV corrections are thus very small here. At high temperatures, 
however, the l/N term is too large for all considered fields. Although the results with the exponentiation trick are 
better and show the correct qualitative behavior, they are not quantitatively better than the MF results. We discuss 
the results further at the end of the next section. 



III. O(N) MODEL 
A. General considerations 



In the last section the Heisenberg model was rewritten in terms of Bose fields and the resulting SU(2) model with 
two boson flavors was generalized to SU(7V). The homomorphism between the groups SU(2) and 0(3) opens another 
way to obtain a large N theory. In this section the Heisenberg model is mapped onto an 0(3) model, which is then 
generalized to O(iV). Many concepts are identical to the SU(7V) case. However, the O(N) model adds a number of 
complications. I— . 

The group theory involved here can be found elsewhere.Ej In brief, the Lie groups SU(2) and 0(3) have the same 
algebra, up to different representations. This means that the infinitesimal generators of SU(2), namely the Pauli 
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matrices, have the same commutation relations as the three infinitesimal generators of 0(3), (Xk)i 



-2ie 



ijk- 



Put 



informally, SU(2) and 0(3) have the same local structure, although the global structure is different. Note that we 
could also talk about the SO (3) model instead of 0(3) since the two have the same algebra. 

The upshot of this is that we can map the Heisenberg model onto an 0(3) model. Wc introduce three Bose fields 
b a and let S = —ickapblfip, k — x,y, z, where we again assume summation over repeated indices. It is easily shown 
that the commutators of the spin components S are correct. To restrict the Hilbert space to the physical states 
two constraints are needed, b^ a b a — S and b^b^ = 0. The second constraint needs explanation. Let us consider the 
eigenstates of S z for a single spin. Since S z is not diagonal in boson flavors we introduce new bosons c\ = (6i +i&2)/v / 2, 
C2 = (bi — i62)/v / 2j and C3 = 63. Then S z — —c\ci+c\c2 and the constraints read c^Cq = S and cldj + c^cj + CgCg = 0. 
The eigenstates of S z are simultaneous eigenstates to cjci, <4c2, and C3C3. As an example, the following table shows 
the eigenvalues of the number operators and of S z for S = 2. The first constraint means that there are two bosons. 



c\ci 


4 C 2 


4 C 3 


s z 


2 








-2 


1 





1 


-1 


1 


1 











1 


1 


1 





2 





2 








2 






The state with S z = is obviously counted twice. The second constraint just removes the last state. It can be 
rewritten as CgCgl^) = — 2c|c2 t where is any state. This means that the state one gets by creating two C3 
bosons is the same as the one produced by creating one c\ and one C2 . The second constraint thus reduces the Hilbert 
space by identifying states with one another. For general S, the second constraint removes spurious spin multiplcts 
of lower total spin. 

The first constraint does not make sense for half integer spin. We assume S integer. We will not have to do this 
for even N in 0(N) theory. 

The 0(3) spin matrix should be an element of the algebra, which consists of antisymmetric 3x3 matrices. In three 
dimensions any antisymmetric matrix is dual to an axial vector. This is in fact the reason why angular momenta 
can be written as axial vectors in three dimensions. Here we go the opposite way and define the spin matrix by 
Sp = ie a/ skS k — b\j)p — b'^b a . Using the antisymmetry of Sp, the Hamiltonian (|l]) becomes 



H 



(73) 



with h = ((0, i, 0), H, 0, 0), (0, 0, 0)). 

To generalize the model to O(N), we introduce N Bose fields b a subject to the constraints 



bib a 



NS 

bibt = 0. 



(74) 
(75) 



The O(N) spin matrices are S9j = b^bp — btb a . The second constraint again restricts the Hilbert space by identifying, 



say, b N b^ N \tfj) with another state. The O(N) Hamiltonian is 



^ = -|^EW5Cj)-|E^«) 



(76) 



where h contains N/3 copies of the 0(3) matrix along the diagonal. NS be an integer multiple of 3. 

The next steps are similar to the SU(A^) model. Going over to the continuum and inserting bosons we get 



H= d 



O T 

js(d 3 bi)(d 3 b a ) - —bi{d 3 b ] p )bp{d 3 b a ) 



B 

-^hfjbpbu 



(77) 
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In writing the partition function as a functional integral, the first constraint ( |74| ) is implemented using a Lagrange 
multiplier field A. Two Lagrange multipliers fi\ and 1x2 are introduced to enforce the two components of the second 
constraint (f75|). They couple to the b a fields in the form /ii Reb a b a + /j.2 lm.b a b a = /z* b a b a /2 + /i&*6*/2, where we 
have introduced fi = fi± + ifi2, which is somewhat misleading, though, since both fj,\ and [12 have to be integrated 
along the imaginary axis. The partition function reads 



\J^dr J d 2 rC[b;X,^ , 



Z = I D 2 b a DXD 2 fie^q>[ / dr / d 2 r C[b\ A, fj] ] , (7S) 

where 



C = ^b* a d b a + JSidjblWjbc) - ^bUdjb^bpidjba) 

B N S 
-hpbpba + \b* a b a —A 

d o 

+ \^b a b a + l -^b* a b* a . (79) 
The quartic term is decoupled using a Hubbard-Stratonovich transformation, 

C! = ^b* a d b a + JSidjbDidjba) + 3NJQjQj 

+ MJQjbXdjba) - MJQjidjbDbv - ^h%bp a 

+ Xb* a b a -^X + ^*b a b a + ^b* a b* a , (80) 

where Qj is real and a gauge field. As compared to SU(iV), additional complications arise since under gauge changes 
fi transforms like a charge 2 particle, as discussed below. We choose the transverse gauge, djQj = 0. 

B. Mean field theory 

Again, MF theory is exact for N — > 00. We make a static assumption for A, Q, and /1. We then express the fields 
b a in terms of Fourier transforms. The partition function reads 

Z = J D 2 b a (k, iuj n ) exp ( - 3MN/3 JQ ■ Qa 2 

+ jtf a **™0X-J d 2 fc]T (81) 



with 



C'a[b] = Pa 2 ^ (-itWn + JSk 2 a 2 - 6 JQ • ka 2 + a 2 A) 6* (k, iu n )b a (k, iu n ) 

a 

a 2 

- Pa 2 2J Blip bp(k, iuj n )b a (k, iui n ) + (3a 2 2J yM* b a (-k, -iui n )b a (k, iui n ) 

a/3 a 

+ /3« 2 E T^" k ' " lw «)^(k, iu n ), (82) 



where sums are again written out. To diagonalize the Zeeman term we substitute new fields, 

C3n+1 = ~ /=(&3n+l + *&3n+2)> = (&3n+l _ *&3n+2), C3n+3 = ^3n+3- (83) 

After shifting the integration variable k to k + 3Q/5 we get 
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4'[c] = /3a 2 (~iHu)n + JSk 2 a 2 - Bh a a + a 2 \ - • QaA c* (k, iw n )c a (k, iw n ) 

a 2 a 2 

+ ^ X! ~2~P* C ^(~ k ' -*^n) c a( k : iw n) + fa 2 ^ ~~^~~P C w(-K -«W n )c* (k, iw n ), 
a a 

where /i is diagonal with the diagonal elements —1,1,0,— 1,1,0,... and 

( 3n + 2 for a = 3n + 1 
a = < 3n + 1 for a = 3n + 2 
3n + 3 for a = 'in + 3. 



(84) 



(85) 



The partition function depends on A and Q only through A = a 2 (3\ — (9/3J/S)Q ■ Qa 2 . To get rid of the terms mixing 
(k, iu n ) with (— k, —iu n ) we note that c a (k, iu n ) is even in w„ and introduce new fields, 



d a (k,iuj n ) = -^=[c a (k,iu; n ) - ic a (-k, iu n )]. 



(86) 



Then we have 



£%[d\ = f3a 2 J2 (~i^n + JSk 2 a 2 - Bh° + ^) d* (k, ioj n )d a (k, ioj n ) 

■ 2 - 2 

+ (3a 2 22 -7T~P* d a(K iu n )d a (k, iuj n ) - /3a 2 2J —]Id^(k, iu n )d* a (k, iu n ). 



(87) 



The fields d a are now integrated out. We define J = (3 J and _B = /3B (note different definition of B). By integrating 
over the d a , putting the product into the exponential, and evaluating the flavor sum we obtain 



Z oc exp 
+ ln 



. JVS T Afa 2 N 

M ~ A - i^y 



In 



{-i(3?iu) n + JSk 2 a 2 + A) - a A f3 2 JI*-p 



(-i(3hu n + JSk 2 a 2 + A) 2 - a 4 (3 2 Jl*-p - B 2 



The MF values A and n are determined by the saddle-point equations 



9A 



ainZo Q ainz _ Q 



The last two are equivalent. They yield 



Q = a^ 2 JI fd 2 k^2(\ 



djl* 



1 



{-i(3hu n + JSk 2 a 2 + A) 2 - a 4 /3 2 7I*7I 
1 



{-i/3hw n + JSk 2 a 2 + A) 2 - a 4 f3 2 JI*-p - B 2 
One solution is 71 = 0. For 71 ^ we get 

1 = [l - 2e^cosh(a 2 /3|7i|) + e~ 2X 



1 - 2e~ A cosh A/a 4 /3 2 |7J| 2 + B 2 + e 



-2A 



(88) 



(89) 



(90) 



(91) 



which has solutions for \a\ > 0, corresponding to broken gauge symmetry. The saddle-point equation for A yields 



l-2e- A cosh(a 2 /3|7i|)+e- 2A 
1 - 2e- A cosh \ja 4 f3 2 \ji\ 2 + B 2 + < 



-2A 



(92) 
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so that there are no simultaneous solutions with ~p ^ 0. Note that if we had ~p ^ 0, a term like \ {dj — 6Qj/S)[j,\ 2 would 
appear in the gauge invariant Lagrangian, which would make the gauge field Q massiveEj (this is the Anderson-Higgs 
mechanism). In our case, however, it is massless at the saddle point. 
The partition function is now 



Zq oc exp 



The MF equation for A become: 
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and the MF magnetization is 



.NS- 



/^EE 



x In ( -iphu n + JSk 2 a 2 - Bh a a + A 



S 



1 



ln(l 



-A+B 



) + ln(l-e- A ) 



4ttJ5 
Hn(l- e -") 



(93) 



(94) 



M = 



MN(3 dB 
1 



lnZ 



ln(l-e- A+a )-ln(l-e 



-A-B\ 



AttJS 

which exhibits the same universality as the SU(iV) result. At low temperatures, 

1 . 1 



M -5 = 



AttJS 



ln(l 



-pB\ 



2irJS 



ln(l 



2/3BN 



(95) 



(96) 



up to exponentially small corrections to the magnetic field. Thus, although the leading term is the same as the 
non-interacting magnon approximation, Eq. ((2^), the second term is different. Thus we expect the correct behavior 
at the lowest temperatures but deviations already for T ~ 2B. 



C. 1/N corrections 

The method used to calculate the 1/N corrections is similar to the SU(A r ) case. However, the second constraint 
([75|) introduces additional problems. The magnetization is now 

m = !e^m = |fE^ c «>' ( 9? ) 

a/3 a 

using the definition of c a in Eq. (|S3|). Fourier transforming the c a we find (c^Cq) = (ct^da) (no summation implied) 
and 

M = |?E^<4<2«}. ( 98 ) 

a 

In the following we use the representation in terms of fields d a , Eq. ( |86| ) . The fluctuations are written as 

A(r, -r) = A + iAA(r, r), (99) 

Hi(t,t) = + iAMi(r,r), (100) 

fi 2 (r,T) = + iA M2 (r,r), (101) 

Q j (r,r) = + AQ j (r, T ), (102) 

where AA, AQj, A/ii, and A/^2 are all real. For convenience we use a complex A/x = A^ii + iA[i2 so that /x(r, r) = 
iAn(r,r) and (i*(t,t) = -iA/i*(r,r). The SXJ(N) methods of Ref. |l7|can be adapted to the O(N) model; we write 
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Z = J DA\D 2 AfiDAQ J cxp(-NS) (103) 

and expand the action S as in Eq. ( |30| ) for SU(iV), where ri can also stand for A/i or A/i*. We can also write 
S = So + Sdii + Sioop with 

So = ^TrlnGo\ (104) 
1 f m f 2 / ~ ~ NS 



Sdir = J dr J d'r \3NJQ • Q - —A ) , (105) 



Sloop = Tr In ^1 + Go ^ 



(106) 



The (normal) MF Green function can be read off from Eq. (|87j), 

G£(k, £w n ) = {-i(3fiuo n + JSk 2 a 2 - Bh a a + X) -1 . (107) 

For /I = no anomalous Green function exists since there is no d* a d* a term. 

The constant part of Sdir together with So again gives the MF action S™ . The first-order terms cancel. The vertex 
factors can be found in analogy to SU(iV), 

^a Qj t = ^(-6J)a 2 %. (108) 

The diagrammatics are similar to the SU(iV) case. For AA and AQ fluctuations the only differences are: (i) The AQj 
vertices contain an additional factor of 3 each, giving 9 in S^ , (ii) the direct AQj propagator from Sdh- contains an 
additional factor of 3, (hi) h% is replaced by h%, (iv) now B = j3B, and (v) A is given by Eq. (E3). 

(2) — 

In particular, we find <SXq 2 (o o) aq 2 (o o) = as for the SU(iV) model. Thus gauge fluctuations are massless for 
O(N) as well as for SU(iV). 

The contribution from Afj, requires some thought. From Eqs. ( |S0| ) and (|8^) we see that A/i couples to two "creation 
operators" d* a d^, whereas A/i* couples to d a d^- Consequently, the boson loop in S^ can only contain one A/i vertex 
and one A/i* vertex or neither of them. Thus <S' 2 ' and the RPA propagator do not mix Afj with other fluctuations. 
Consequently, the only contributions to (c£t <f Q ) at the 1/N level correspond to the diagrams in Fig. [?], where the 
zig-zag line denotes the Afj, RPA propagator. Note the directions of the boson lines. 

We derive sj^, Ap , S AX , Q ^ A r A , and S^t^, A ^ in App. ||. We can then integrate out the fluctuations. Afi is 
a complex field so that the contraction of a pair yields 



z f D 2 Zt zlz £2 e^- Y S^,z e , j = -(5^)^. (109) 
Consequently, the diagrams of Figs, ||(a) and 0(a) added together, and Figs. ||(b) plus 0(b), respectively, are 

(110) 




„ o jvn/6j + j: p ^ Epd, 



n B {K) + n B (K + B) 

°- (Ml. 



q ,^ \EX ' m/Qj + T.poi EX 



where 




n B( e £+q/2) + n s( e £- q /2) 



(112) 



2JSk 2 a 2 + JSq 2 a 2 /2 + 2A 
4 ) (^ +q/2 )+4 1) (eg_ q/2 ) 2 l + "B(^ +q/2 )+"B(^- q /2) \ f . 



-i(3hv n + 2JS k 2 a 2 + JS q 2 a 2 /2 + 2A {-i(3hv n + 2JS k 2 a 2 + JS q 2 a 2 /2 + 2A) 2 
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the other symbols are identical to the SU(iV) case if h% is replaced by h™. 
With Eq. ( p8|) the 1/N contribution to the magnetization reads 



■ ' ■/ trv ex nn/6j +j2/3 a ± ex / 



TV 8vr 2 
with 



n B (A-B)-n B (A + B) 

ci = = = = = — . (115) 

n B (A - B) + ns(A) + n B (A + B) 

Evaluation of the AA and AQ2 contributions is analogous to the SU(Af) case. In particular, naive summation over 
iv n results in a strong divergence. The constant term in the integrand for large momenta is 3(ci — Mo/ S). Numerical 
calculations confirm this result. Again, the spuriou s div ergence is removed by taking operator ordering into account. 

We now turn to the A/j, contribution. From Eq. (112) we see that cr" diverges logarithmically at large momentum 
k because of the summand 1 in the numerator. However, the A^i contribution to the magnetization is finite. To see 
this, we use a finite cutoff K and let K — > 00 in the result, cr" is dominated by 



ol = 



f d 2 k 



4tt 2 Jk<K -i(3hv n + 2JSk 2 a 2 + JSq 2 a 2 /2 + 2A 
A hn{-if3hv n + 2JS K 2 a 2 + JS q^a 2 /2 + 2A) - ln(-i{3Kv n + JSq 2 a 2 /2 + 2A) . (116) 



8ttJS 

For a"' the corresponding contribution is 



8irJS V -iPKv n + 2JS K 2 a 2 + JS q 2 a 2 /2 + 2A -i/3hv n + JS q 2 a 2 /2 + 2A, 
Note that these two expressions do not depend on a. The frequency sum over the A/i contribution is 

1 1 



(117) 



1 ]n(-i/3hv n + e K ) - ln(-i0nv n + e ) ' 

where tk = 2JS k 2 a 2 + JS q 2 a 2 /2 + 2A. Since this contribution is proportional to c\ it comes only from the diagram 
Fig. 0(b). 

The sum over iv n can be evaluated by contour integration. As noted in App. 0, splitting the time to enforce 
correct operator ordering results in an overall factor of exp(iv n ri) , which removes any ambiguity in the iv n sum. In 
the complex v plane, AM* has a branch cut along the real axis between the points eo/h[3 and e^/hp and two poles on 
top of the branch points. The contour integral contains three terms: Two from integrating around the branch points 
in small semicircles and one from integrating the spectral function of AAf* along the branch cut. The two semicircles 
contribute — 2ci[ns(eo) + ub^k)]- For K — + 00 this expression becomes — 2cins(eo). The spectral function is 

P- 1 



Anci ex/hp— v is-e /hj3 

W \\n(e K /h0-u)-hi(v- co/hptf 



A = _z^i_ — - zz 2 ( 119 ) 



and the integral over it can be shown to vanish for K — * 00. Thus 

AM* = -2cm B (j Sq 2 a 2 /2 + 2A) (120) 



for K — > 00. With Eq. (114) the full contribution from A/i to the magnetization is 
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n B (K-B) -n B (K+ B) 
n B (A -B) + n B (A) + n B (A + B) ' 



(121) 



A few remarks are in order: (i) This contribution increases the magnetization, whereas fluctuations in A and Q 
decrease it. The physical explanation is that the MF approximation, which enforces the second constraint b^b^ = 
only on average, underestimates the magnetization because it contains contributions from spurious multiplets of lower 
total spin, (ii) The A^i contribution has a typical energy scale of 2A since excitations of energy 2A (and higher) are 



removed by the second constraint (75). (iii) The q integral over AM* is well-behaved for large q so that a cutoff, 
which is necessary for the AA and AQ contributions, does not change the result appreciably but would complicate 
the calculations. 

Figure H shows the magnetization for S — 1/2 and B/J = 0.05,0.1,0.25 as a function oiT/B. The A/j, fluctuations 
win over the other contributions; the magnetization is larger than the MF result. We see that the 0(A) 1/N expansion 
gives much better results than the SU(A) model except at low temperatures. At small T/B the magnetization seems 
to be unphysically large, especially for smaller fields. At moderate temperatures the 0(A) 1/A magnetization is 
better than both MF results. At high temperatures, the Monte Carlo data consistently fall slightly below the O(A) 
1/N and, for B/J = 0.25, even below the O(A) MF results. 

Recall that the SU(A) MF approximation works well at low temperatures because it coincides with the non- 
interacting magnon approximation for the Heisenberg model, whereas the 0(A) MF magnetization does not. The 
1/A corrections for the 0(A) model (predominantly from A/x at low T) are large and in fact overcompensate for the 
error made at the MF level. 

There is a distinct crossover to the moderate T regime, where SU(A) MF becomes too large, SU(A) 1/A becomes 
quite wrong, and 0(A) 1/A is rather good. In fact it is surprisingly good considering that 1/A is not really small. It 
is not fully clear why the 0(A) model works better than the SU(A) model at moderate and high temperatures. The 
reason may lie in the different behavior of gauge (AQ) fluctuations in the two models.Ej In the SU(A) model they are 
massless in general, whereas for 0(A) they are massless only because the MF value of /i happens to vanish for this 
particular system. In both cases the zero mass leads to an overestimate of fluctuations at the 1/A level. However, in 
the 0(A) model fluctuations in /i are available to compensate for this, thereby partly restoring the effect a massive 
gauge field would have. 

One might think that 0(A) should.-be worse since the 0(A) model for A > 3 does not have skyrmions, whereas 
the SU(A) model has them for all A. El However, the 1/A expansion does not contain these iiqn-perturbative effects 
anyway. On the other hand, they are, in principle, captured by the Monte Carlo simulations.113 

The deviations between 0(A) results and Monte Carlo data at high temperatures and B/J = 0.25 or larger (see 
Fig. ^) are probably due to thermally created skyrmions or to the fact that the simulations are done on a lattice, 
whereas the 1/A calculations use a continuum approximation. The dispersion of the former is a cosine band if 
bosonized, whereas the latter has parabolic dispersion. Both effects should become important for temperatures T > J 
since both the band width and the typical skyrmion energy are of the order of J. Indeed, the deviations start at T ~ J. 
(In the same region higher order gradient terms not included in the Heisenberg model should become important.) 

We have also investigated the universal dependence on JS 2 and B. Whereas the MF results-exhibit this universality, 
it is violated by small logarithmic corrections at the 1/A level for both models, as expected.liil 

Our results can be compared with the microscopic approach of Kasner and MacDonald,Lj which includes spin- wave 
corrections to the electronic self-energy. This approach is microscopically better justified than the Heisenberg model. 
However, the magnetization from Ref. [l4]is consistently too large and even MF SU(A) and 0(A) results agree better 
with Monte Carlo data. n 

Comparison to NMR experiments by Barrett et alB shows a number of discrepancies. At low temperature, the 
experimental data look flat, whereas at high T they drop well below the theoretical results. These discrepancies are 
mainly due to the scaling of the data, which is done by setting the measured magnetization, which is reduced by 
disorder, to S in the limit T — ► 0. n 

Recent magnetoabsorption measurements by Manfra et ala show better agreement with our results. In Fig. |^ we 
compare data of Ref. || with SU(A), 0(A), and quantum Monte Carlo results. In the calculations we have used the 
exchange constant corrected for finite width of the quantum well, which yields B/ J w 0.32. The Monte Carlo data fall 
below the 0(A) 1/A results above T ~ J, as discussed above. The experimental data agree quite well with SU(A) 
theory at low temperatures and with 0(A) 1/A (and Monte Carlo) results at moderate temperatures, as expected. 
At higher T the experimental data show more noise but lie mostly above the 0(A) 1/A curve. This discrepancy 
for T > J is probably due to neglected higher gradient terms in Eq. (jl]). The experimental system is a continuous 
itinerant magnet, which probably explains the deviations from Monte Carlo lattice simulations. 
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IV. SUMMARY AND CONCLUSIONS 



We have calculated 1/N corrections to large N Schwinger boson mean field theories for the two-dimensional fer- 
romagnetic Heisenberg model, meant to describe a quantum Hall system at filling factor v = 1. Normal ordering of 
operators has to be carefully taken into account to obtain the corrections. Using a O(N) mode^-iwe find reasonable 
agreement of the 1/N corrected magnetization with both quantum Monte Carlo simulationstSEJ and experiments^ 
at moderate and higher temperatures. At low temperatures, the SU(AT) model works better since it reproduces the 
correct low-energy physics. However, the SU(iV) model does not describe the data anywhere .else, confirming Auer- 
bach's remark that large N methods are "either surprisingly successful or completely wrong. "113 Effects of thermally 
created skyrmions, which are not included in our approach, are small. Away from filling factor v = 1, skyrmions are 
present in the ground state and should be important. The natural next step leading on from this work would be to 
incorporate these skyrmions. In addition higher derivative terms due to the long range Coulomb interaction should 
be investigated. 
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APPENDIX A: CALCULATION OF DIAGRAMS FOR SU(iV) 



Here, we derive explicit expressions for S^ 2 \ S^, S^ 1+1 \ and 5^ 2+1 - ) . We first calculate the Gaussian part 
S^J 2 ri 1 ri 2 /2 in the action. It consists of two contributions, see Figs. [l](b) and ^(b). The first term is easily read off 
from Sdir in Eq. (||), 



1 o(2) 

2 ^2 



dir 



4^ 



^ J AQ(-q, -iv n ) ■ AQ(q,w„ 



(Al) 



Note that AQ 3 (— q, — iv n ) — AQ*(q, iv n ) since AQj(r, r) is real. The same holds for AA, below. 

The loop part follows from Eq. (|3q). The notation is shown in Fig. |Io|. Inserting Eqs. (j34|), (|37|), and (^8|) we obtaii 



c(2) 
2 ^iik 



loop 



r ^ = E E E 

l 



2?r r 



-i/3huj n - ij3hv n + JS(k + q/2) 2 a 2 - Bh% + A M 



i/3a 2 AA(q, iv n ) - 2Ja 2 (k + q/2) • AQ(q, iv n ) 



1 



2tt 



i/3a 2 AA(-q, -iv n ) - 2Ja 2 (k - q/2) • AQ( q, -iu n ) 



-if3huj n + JS(k - q/2) 2 a 2 - Bh% + AJV 
Performing the Matsubara sum over iui n and utilizing the periodicity of the Bose function ns we get 
1 



1 q( 2 ) 
2 



loop 



1 47T 2 ^ "■B( £ k+q/2) - "a( e k-q/ 2 ) 

rt x rt 2 - + 2N N 2 ^ -ij3Hv n + 2 JS q • ka 2 

i/3a 2 AA(-q, -iu n ) - 2Ja 2 (k - q/2) • AQ(-q, -iu n ) 
i/3a 2 AA(q, iv n ) - 2 Ja 2 (k + q/2) • AQ(q, iv n 



(A2) 



(A3) 



where e£ = JS k 2 a 2 — Bh^ + A. At this point we use the transverse gauge, q • AQ(q, iv n ) = 0. We choose coordinates 
in such a way that k\ and AQi are parallel to q. Then A<5i(q, w n ) — and the last expression simplifies, 
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I C( 2 ) 

2 



loop 



re ^ 



1 4tt 2 n B (eg +q/2 ) - n B (eg_ q/2 ) 

2iV ^ v „ -<^n + 2 JS gfcia 2 



/3 2 a 4 AA(-q, -w„)AA(q, w„) + 4J 2 a 4 fc 2 AQ 2 (-q, -w„) AQ 2 (q, w„) 



Terms mixing AA and AQ vanish since their coefficient is odd in fc 2 . Adding Eq. (Al) to Eq. (A4) yields 



4^ 2 



? (2) 

, AA(-q,-w„),AA(q,w„) _A/"2 



EE 



"fl( e k+q/2) - n B( £ £-q/2^ 



-if3%v n + 2JS qkia 2 



5 



(2) 

AQ 2 (-q,-if7 l ),AQ 2 (q,ji'n 



87T 2 



Ja 2 



in 2 



-iftTiVn + 2 J S qkia 2 



2- 



(A4) 

(A5) 
(A6) 



all other components vanish. The fact that only connects fluctuations at (q, iv n ) and (— q, — iv n ") just means 
that the RPA propagator conserves energy and momentum. The real part of S^ 1 is always positive except for 

(2) 

"AQ 2 (Q o) aq 2 (o o) = 0- Thus there is one zero mode, which results in an additional factor in the partition function, 
which, however, does not depend on field and is thus irrelevant for the magnetization. The zero mode at q = 0, 
iv n = shows that gauge (AQ) fluctuations are massless. For the remaining modes, can be inverted to get the 
RPA propagator D, which is also positive. The saddle-point is thus stable. 

Looking at the diagrams in Fig. ^ we see that the horizontal propagator in the right diagram can only be at q = 0, 
iv n = since the source j a does not insert any frequency or momentum. In fact it can only be AA(0, 0), as we will see. 

Keeping this is mind we calculate 5^ 3 ] ^ . Since is symmetric in its indices we can assume that ri 1 is AA(0, 0). 
Furthermore, rg 2 determines rg 3 . We start from the definition (p5|), 



S. 



(3) 



1 

N 



[Tr {GoVt x GoVi 2 GoVi 3 ) + Tr (Gove^ov^Gove^) 



(A7) 



The first of the two summands is 



-i/3hu) n + JS(k- q/2) 2 o 2 - Bh% + A 
1 



-i(3huj n - i/3hv n + JS (k + q/2) 2 a 2 - Bh% + A -if3huj n + JS (k - q/2) 2 a 2 - Bh<* + A 



k 

d 



n B {JS (k + q/2) 2 a 2 - Bh% + A) 
{i[3Tiv n - 2JSqkia 2 ) 2 

n B (z) 



dz - z - i/3hv n + JS(k + q/2) 2 a 2 - Bh% + A 



-JS(k-q/2) 2 a 2 -Bh%+A 



n B {JS (k + q/2) 2 a 2 - Bh a a + A) 
{-il3hv n + 2JSqk 1 a 2 ) 2 



n B (JS (k - q/2) 2 a 2 - Bh a a + A) n { B > (JS (k - q/2) 2 a 2 - Bh£ + A) 



(—i(3hv n + 2JS qk±a 2 



—i(3Tiv n + 2JS qk\a 2 



(A8) 



where n^ B \e) = d v n B (e)/de 1 ' is the ^-th derivative of the Bose function. With the vertex factors the last expression 
becomes 



(27T) 3 

Af 3 N 



EE** 



f -(3 2 a 4 
\ Apa^k 2 



n B (JS(k + q/2) 2 a 2 - Bh a a + A) 



(—ififiVn + 2JS qk\a 2 



n B (JS(k - q/2) 2 a 2 - Bh a a + A) 



(—if3Tiv n + 2JS qkia 2 ) 



2 \2 



(JS(k - q/2) 2 a 2 - Bh% + A) 
—i(3fiv n + 2JS qk\a 2 



where the upper [lower] term in the curly brackets is for rg 2 = AA(q,iz/ n ) [AQ 2 (q, iu n )]. For 
integrand would be odd in fc 2 so that this contribution vanishes. 



(A9) 

AQ 2 (0,0) the 
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The second term in Eq. (A7) just has rg 2 and rg 3 exchanged, which means q and iv n have opposite sign. Thus, 

5 (3) - ( ^ )3 i3d 2 V V ( £ k+q/2) ~ ( e k-q/ 2 ) , 2 4 ^ / A10l 

C, AA(0,0),AQ 2 (-q,- ! i'„),AQ 2 (q^„)" A/"3 /V P Z^ Z^ „-„*„ , 7C„7, „2 2 ' 1 ^ 



where the fractions in Eq. (|A9|) containing the denominator squared have cancelled upon adding the two terms. 
We can now calculate and S^ 2+1 \ These expressions contain a vertex Vj a = 2ita 2 /AT instead of v\\ — 

2iria 2 /Af. The source j a inserts zero momentum and frequency. For <S( 1+1 ) we are thus only interested in Sj 1+ &x( j 

[the left loop in Fig. Q(b)]. By taking the limit to zero frequency and momentum, we obtain 

? (2) - ^Y J Y. n B ) i JSk2a2 - Bh a a +m-p 2 a A ) 



J AA(0,0),AA(0,0) 



k a 



^j-B 2 a A [n B {h-B)+n B {k + B)\. (A12) 



27V JS* 

Keeping in mind that j a couples only to the boson of flavor a we find similarly 

CS(o,o) = -j^p 2 a*n B (A-Bh°). (A13) 

From we infer S^ 2+1 \ 

c(2+i) - (2 ^ )3 3a 2 V "^ )(e k+q/2) ~ n 5 ) ( e k- q / 2 ) 2 4 

<^;AA(- q ,-^),AA (q ,^) - M3N P* Z. _ mVn + 2 J Sqha 2 ^ a >> 
5 (2+l) ~ (2 ^ )3 3a 2 V n ^ )(e k+q/2) - n B ] K-^2) j 2 4 2 



APPENDIX B: CALCULATION OF DIAGRAMS FOR O(iV) 
We start with S^ 2 \ using the notation shown in Fig. [H} In analogy to Eq. (|35|), 



f 



^(q^^AKq,^) = Tr (G0V A ^G* V Afl ) . (Bl) 

Here, one of the Greens functions is the complex conjugate since the line is traversed against the direction of the 
boson propagator. Momentum and frequency of Gq are measured counter-clockwise. We find 

c(2) 3 2 a 4 V V - 

A^^iu^AK^) Af 2 N Z_. _ i(3huJn _ + j 5(k + q/2 )2 a 2 - B/tg + A 

1 



-i/3?ia>„ - J5(k - q/2) 2 a 2 + BhZ - A 



tt 2 2j4 yy 1 + n B( £ £+ q /2)+ n5 ( £ k-g/2) 
M 2 N l -i3hv n + 2JS k 2 a 2 + JS q 2 a 2 /2 + 2A 



(B2) 



with e£ = JSk 2 a 2 - Bh^ + A. Here, we have used that h% + h£ = and the identity n B (-e) = -n B (e) - 1. The 
real part is positive so that the functional integral is well-defined. This kind of expression is known from the theory 
of scattering processes. 

5 55(o,o),A/»»,A /1 can be derived similarly, 
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s 



(3) _ 

AA(0,0),A/i*(q,ii>„),A/i(q,ii/„) _A/3jY 



(2tt) 3 i/3 3 a 6 



EE 



? (1) r f tt i 



l B l e k-q/2^ 



k m v -i/3hv n + 2JS k 2 a 2 + JS q 2 a 2 /2 + 2A 
- n B( £ k+q/ 2 )+ n B(e£L q/2 ) \ 



(- j/37iz/ n + 2 JS* Pa 2 + JS" q 2 a 2 /2 + 2A) 2 



(B3) 



The term containing the denominator squared does not cancel in this case. To obtain this result we have summed 
over boson frequencies iuo n using contour integration. One has to consider operator ordering to do this properly. Since 
the anomalous combinations d^d^ a and dad a in the Hamiltonian contain two commuting operators time splitting is 
not necessary and no phase factors appear in and ua m . . On the other hand, the other vertices obtain factors 



exp(iw„7/). It can be shown that the two terms in coming from the symmetrization in Eq. (A7), obtain factors 
of ex jp(im n r)) exp(iv n ri) and exp(— iu> n rj) e~xp{iv n rf) , respectively. The different factors in iui n are crucial in arriving at 
Eq. ( |B3| ). Furthermore, we obtain an overall factor of exp(iiy n rj). 
Immediately we find 



S 



(2+1) 



(27r) a (3 3 a 6 
- 2- 



ja ; A/i* (q,iv„ ) ,A^(q,iv n ) 



E 



n 



S ) ( £ k+q/2) + n S ; ( 6 k-q/2 



(1)/ 



k v -i(3nv n + 2JSk 2 a 2 + JS q 2 a 2 /2 + 2A 
1 + "s(^ +q/2 ) + n B (6«_ q/2 ) \ 



{-i(3hv n + 2JS k 2 a 2 + JS q 2 a 2 /2 + 2Af 



(B4) 



Finally, we have to recalculate S 



,(2) 



(2) 

AA(0,0),AA(0,0)' 



^AA(0,0),AA(0,0) 



By replacing h% by h% in Eq. ( A12 ) we get 



'AAfJS 



^/Ta 4 n B (A- B)+n B (A)+n B (A + B) 



(B5) 
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(a) (b) 

FIG. 1. Diagrams contributing to <Sdi r 




FIG. 3. The diagram for the MF magnetization. 
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(a) (b) 

FIG. 4. \/N diagrams contributing to the magnetization. Diagram (a) corresponds to T>^ in Eq. ( |45| ) and (b) corresponds 
to V { ° ] in Eq. @. 
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FIG. 5. Contour of integration for the Matsubara sum over iv, 
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FIG. 6. SU(iV) magnetization for magnetic fields (a) B/J = 0.05, (b) B/J — 0.1, and (c) B/J = 0.25 as a function of 
temperature in units of B. The thin solid line is the MF magnetization, the thick solid line includes 1/N corrections linearly, 
the thick dashed line includes themiis the exponential to force positivity, and the circles with error bars are quantum Monte 
Carlo results for a 32 x 32 latticeJBEl 
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FIG. 8. O(N) magnetization for magnetic fields (a) B/J = 0.05, (b) B/J = 0.1, and (c) B/J = 0.25. The thin solid line 
is the MF magnetization, theiithick solid curve includes 1/N corrections, the circles with error bars are quantum Monte Carlo 
results for a 32 x 32 lattice|i2Hi3 and the dashed curve shows the SU(A r ) MF magnetization for comparison. 




FIG. 9. Comparison of SU(iV) aiwl O(N) results, as well as quantum Monte Carlo results for a 16 x 16 lattice, with 
experimental data from Manfra et aln. The open squares were obtained sweeping the temperature at fixed field and the filled 
squares by sweeping the field at fixed temperature. 
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k- q/2,ico„ 




FIG. 10. Notation of momenta and frequencies for with external AA or AQ legs. 
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k+q/2,ico„+iv„ 

FIG. 11. Notation of momenta and frequencies for .S*- 2 ' with external A/i legs. The momenta and frequencies are measured 
in the direction of the attached arrows. 
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